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1 . 0 INTRODUCTION 


1.1 Statement of Problem 

One atmospheric phenomenon that may adversely affect the operation 
of aeronautical systems is the cold surge or gust front accompanying 
thunderstorms. As the cold downdraft associated with the mature or 
decaying thunderstorm cell reaches the ground and spreads out horizontally, 
this cool, relatively dense air pushes under the warmer, lighter ambient 
air. The interface between the two air masses of different densities is 
called the pseudo cold front or gust front. This front may precede the 
main storm cell by as much as ten miles. The strong vertical and horizontal 
wind shears associated with the gusty winds behind the front are potential 
dangers to the safe operation of aircraft and launch vehicles. 

In this report we investigate the structural properties of six gust 
fronts. Of particular interest is the variation of wind speed with time 
and height. In an effort to develop a statistical prediction scheme for 
the wind speed from the surface to about 150 m during the passage of a 
gust front, we have scaled the winds so that certain similarities among 
storms are present. We then attempt to relate the scaling parameters to 
measurable synoptic scale and radar variables, based on data from 81 
thunderstorm cases in Tampa, Florida. 

* 

> 

1.2 Previous Research 

The low-level variation of wind under high wind conditions has been 
the subject of numerous empirical studies. As early as 1937 > Sherlock and 



Stout (1937) investigated the effect of wind loading of electric power 
lines during high wind situations. They analyzed the vertical variation 
(up to 250 feet) of wind during two winter storms, and noted that gust 
maxima occurred in the upper levels first: a time lag of several seconds 
existed between the top and the 50 foot level. 

Considerable attention has been directed toward the prediction of 
peak surface wind gusts during thunderstorms. Fawbush and Miller (1954) 
developed an empirical technique for forecasting maximum surface wind speed. 
Their method, called the downrush temperature technique, was based on 
Brancato's (1942) discovery that temperatures of downdrafts reaching the 
surface in thunderstorms are approximately equal to the surface tempera- 
ture of the moist adiabat through the intersection of the environmental 
wet-bulb curve and the zero degree isotherm. 

Foster (1958) devised a method based on the buoyancy equation to 
compute the vertical velocity in a downdraft. This calculated downdraft 
velocity was then correlated with peak wind gusts at the surface. The 
correlation coefficient of 0.5 was disappointingly small, but it was 
statistically significant at the one percent level. To test the 
hypothesis that the downdraft transported upper-level momentum to the 
surface, the average wind speeds at 500 and 700 millibars were added to 
the computed downdraft velocity. This addition yielded an improvement 
in the correlation coefficient of only 0.01. Therefore, little useful 
information had been added. 

Feteris (1965) extended Foster's and Fawbush and Miller's techniques 
by considering combinations of parameters rather than single parameters 
as predictors. An important parameter was the instability index, AT, 
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defined as the difference between Fawbush and Miller’s downrush temperature 
and the ambient air temperature at the surface. When AT was less than 
zero, a strong relation existed between maximum wind speed at the surface 
and mean vertical wind shear in the troposphere. This relationship 
supported the hypothesis that transfer of momentum from the. middle and 
upper troposphere occurred under convectively unstable conditions. The 
results of the multiple linear regression analysis indicated that the 
combination of AT and the surface geostrophic wind speed explained the 
largest percentage of the variance. On the other hand, the predictors 
suggested by Fawbush and Miller explained very little variance of the 
peak wind gust. 

Miller (1967) stated that the downrush temperature technique was not 
designed for tropical-type thunderstorms. Because the wet-bulb zero degree 
isotherm intersection in the subtropics during the summer is well above 
10,500 feet, air from this level seldom reaches the surface. Consequently, 
this technique gives forecasts of wind speed which are too high. Miller 
also indicated that the dry instability index showed promise in predicting 
maximum gust speeds in regions where thunderstorm development is primarily 
associated with surface heating. 

In 1965, data from a micrometeorological network of 14 towers 
equipped with wind sensors at various levels were used to test the down- 
rush temperature technique at Cape Kennedy, Florida (Aerospace Review, 
1965). Although this technique tended to overpredict peak wind gusts, 
it did prove to be useful as a forecast aid. The results from this study 
also pointed to the complexity of the problem of predicting peak winds 
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associated with individual thunderstorms. The regions of high winds were 
extremely localized and sharply edged. Thus a very dense observational 
network is necessary in order to accurately verify prediction methods. 

Hiser, Senn and Courtright (1970) attempted to relate strong winds 
in tropical thunderstorms to radar parameters. Peak gusts were correlated 
with echo features including cell movement, age of storm, horizontal and 
vertical size, growth rates of the thunderstorm top and intensity of radar 
return. Although no single echo parameter could be used to identify 
all severe wind cases in south Florida, a combination of parameters showed 
promise. 

Laboratory experiments have provided a useful analogy between 
atmospheric gust fronts and density currents. Simpson (1969) found 
similarities between the intrusion of a saline solution into fresh water 
and thunderstorm outflows. In most gust fronts examined, internal Froude 
numbers (ratios of inertial to buoyancy force), were similar to those found 
in tank experiments. 

Colmer (1971) summarized the characteristics of gust fronts in 
different parts of the world, including several examples from Oklahoma. 

The total wind vector change associated with Oklahoma gust fronts was much 
greater than in the Bedford, England or Belgian Congo gust fronts, due 
to both the greater speed and greater direction changes. Analyses of the 
vertical cross sections of wind speeds of 13 gust fronts indicated that 
only five resembled Simpson's laboratory density current model which 
showed a lag near the lower boundary. Colmer suggested that this lag occurs 
only intermittently, being destroyed by cold dense air from aloft, then 
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reforming due to surface stress. The average gust size (wind speed 
increase) showed an increase with height in the lower layer (about 500 m) , 
with most of this increase occurring in the first 50 m. The average 
gust length (the distance over which the gust size occurred) showed a 
slight decrease with height. Wind direction shifts were nearly constant 
with height in most of the 13 gust fronts. 

Larger gusts were associated with faster moving storms. Differences 
in the gust fronts were also related to differences in the age of the 
storm. The gust front tends to occur close to the storm center while 
the storm is in the developing and mature stages, and then moves further 
ahead during the decaying stages. On the average, the wind shift line 
preceded the gust front by only a few minutes. The gust size varied in 
a direction parallel to the gust front, with the maximum occurring close 
to the projected storm track. 

A three-dimensional analysis of a severe gust front associated with 
the Oklahoma squal line of May 31, 1969, was presented by Charba (1972). 

A windshift and pressure jump preceded the cold air outflow. Charba 
suggested that this pressure discontinuity was the result of a gravity 
wave jump. A projecting nose, similar to those found in laboratory 
currents, occurred about 750 m above the ground. It was also demonstrated 
that the gust front was dynamically similar to theoretical and experimental 
gravity currents by the near-equal values of the internal Froude numbers. 

The results of previous investigations that have attempted to 
develop prediction schemes for the peak wind gust in a thunderstorm have 
been only partially successful in isolating those processes or parameters 
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that are important in the determination of the gust speed. Even the best 
results explained only sixty-six percent of the variance (Feteris, 1965). 

This disappointingly low value points to the complexity of the thunder- 
storm and the associated thunderstorm gust front. The circulations associated 
with thunderstorms are fully three-dimensional, non-linear, and non-steady. 
Energy sources and sinks are large. The primary causes of thunderstorms 
may differ from storm to storm. The environments may also be quite 
different. Therefore, it is not suprising that no single parameter explains 
very much of the variance in thunderstorm peak wind speeds. 

1.3 Specific Goal of This Research 

In order to predict accurately the maximum wind speed associated 
with any given thunderstorm, it would probably be necessary to develop a 
time-dependent, three-dimensional model capable of interacting with its 
environment. Such a model would require a representation of the important 
cloud physics as well as the cloud dynamics. Development of such a model 
would require an extremely large expenditure of money and effort, and is 
clearly outside the scope of this thesis. Instead, we investigate the time- 
dependent vertical structure of three Florida gust fronts and compare these 
subtropical storms with three continental severe storms in an effort to 
ascertain whether any simple, general properties of the wind structure exist. 
We further attempt to relate statistically any general properties of the 
wind structure to measurable synoptic and/or radar parameters utilizing 
data from Tampa, Florida. 
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2.0 THE OBSERVATIONS 


2.1 Kennedy Tower Data 

NASA’s 150 m meteorological tower at Kennedy Space Center is 
located within Complex 39 Area on Merrit Island. There are no man- 
made structures or natural obstacles within several hundred feet of 
the tower which would cause unrepresentative distrubances in the ambient 
air flow. Wind speeds and directions were measured at 0.1 second 
intervals by means of dual-mounted sensors with the aid of a wind 
selector which automatically selected the proper side of the tower for 
recording. The data were collected at heights of 18, 30, 60, 90, 120 and 
150 m and recorded on magnetic tape. Additional information concerning 
the 150 m tower is given in NASA Technical Memorandum X-53699. 

The recorders were activated only when it seemed likely that a 
thunderstorm would pass over the tower. From the recordings, one- 
second averages of the data were computed and time-height cross sections 
were constructed. An example of a one-second average, time-height cross 
section is shown in Figure 1. Wind information was made available for 
21 situations. Of these, only three cases had strong, well-defined 
maxima which could be positively identified with a thunderstorm cell on 
the radar film obtained from Daytona Beach. Due to both the lack of 
temperature data and the great distance (50 mi) between Cape Kennedy and 
the Daytona Beach radar, it is impossible to say definitely that any or 
all of the three Kennedy wind increases are associated with gust fronts. 
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Figure 1. Time-height cross section of the one-second 
average wind speeds for the Florida gust 
front of July 27, 1967. 
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However, in this study, all wind speed increases such as those contained 
in these cases shall be referred to as the speed maxima associated with 
gust fronts. 

2.2 NSSL-WKY Tower Data 

To compare the gust fronts associated with Florida thunderstorms 
with more severe midwest thunderstorms, three severe thunderstorm gust 
front wind records were obtained from the National Severe Storms Laboratory 
at Norman, Oklahoma. Wind data for these three cases were recorded at 
the WKY television tower which is also used by the National Severe 
Storms Laboratory (NSSL) for meteorological data collection. This tower 
is located approximately seven miles north-northeast of downtown Oklahoma 
City and 23 miles north of Norman, Oklahoma. The surrounding area is 
relatively undeveloped. The wind sensors are Bendix-Friez aerovane 
wind transmitters mounted on booms at heights of 7, 44.5, 90.5, 177.1, 
266.2, 355.4, and 444.6 m. The signals from the aerovanes are recorded 
by analog wind recorders operating at six inches per hour. A more detailed 
discussion of the WKY tower is contained in the NSSL Technical Memorandum 
#50. 


2 . 3 Tampa Data 

The three Florida cases and the three Oklahoma cases comprise the 
gust front sample from which the vertical and temporal variation of the 
wind speeds during the passage of a typical gust front is inferred. It 
will be shown that wind information near the surface is sufficient to infer 
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winds at higher levels. 

The second part of this report is concerned with relating parameters 
near the surface which describe the gust fronts, such as the maximum wind 
speed, to synoptic scale and/or radar parameters. No tower is required 
for this analysis. Data were obtained from the National Climatic Center 
and were used in statistical analysis to determine the predictability 
of these descriptive parameters. These data consisted of surface weather 
observations, radar observations, surface wind recorder charts (1200Z 
and OOOOZ of the following date) , winds aloft and rawinsonde temperature 
and humidity information (1200Z and OOOOZ of the following date) for 81 
storms at Tampa, Florida. The data spanned the summers (June through 
September) of 1967 through 1971. 
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3.0 ANALYSIS OF TOWER DATA 


3.1 Relevant Time Scales 

Before describing the analysis of the tower data, it is important to 
define the relevant time scales associated with the gust front. In general, 
atmospheric circulation features are categorized by three main time scales, 
the synoptic scale of a few days, the mesoscale of a few hours, and the 
microscale of a few seconds to minutes. From inspection of many times 
series of wind speeds during the passage of a thunderstorm gust front, 
there appear to be at least two distinct time scales. One is associated 
with the microscale turbulence and has a charateristic period of a few 
seconds. The other is associated with the organized longer-lived circulation 
of the thunderstorm itself and has a typical period of about 20-30 
minutes. This period is obviously related to the time scale associated 
with individual convective elements. 

The latter time scale is of primary interest in this study, For 
example. Figures 2a, 3a, and 4a show 130-second average time-height 
cross sections of wind speed at Cape Kennedy during the passage of three 
thunderstorm gust fronts. The prominent feature in each figure is the 
wind speed jump from a minimum to a maximum. Associated with this increase, 
AV, is a time interval. At. The average At for the wind ijicrease in all 
six gust front cases was about 11 minutes. The values of At and AV are 
believed to be the most important characteristics needed for gust 
prediction. The higher frequency oscillations will be considered essent- 
ially random and unpredictable. 
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Figure 2a. Time-height cross section of the 130-second running mean 

of wind speeds for the Florida gust front of July 27, 1967 
The isotachs in this figure are expressed in mps. 
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Figure 2b. Time-height cross section of the 130-second running mean 

of wind directions for the Florida gust front of July 27, 1967 
The isogons in this figure are expressed in degrees. 
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Figure 4b. Time-height cross section of the 130-second running mean of wind directions 
for the Florida gust front of June 25, 1968. The isogons in this figure 

are expressed in degrees. 



3.2 Smoothing the Wind Data 


The original wind data for the Florida cases consisted of data 
at 0.1-second intervals stored on magnetic tape, which had been 
averaged over one-second periods. To isolate the variations in wind 
associated with the longer time scale, 130-second moving averages 
were computed every ten seconds. 

The original Oklahoma data consisted of strip charts of wind traces 
recorded at the rate of six inches per hour. These traces were first 
smoothed manually. Time-height cross sections were then constructed by 
selecting wind speeds at 2.5-minute intervals. These analyses are 
presneted in Figures 5, 6, and 7. 

Although a sampling interval of 2.5 minutes is rather long to 
represent a phenomenon with a period of about 20 minutes, prior smoothing 
of the time series should reduce the effect of aliasing on the analyses. 

The smoothed time-height cross sections (Figures 2-7) formed the 
basis for subsequent analysis. A brief qualitative description of the 
six gust fronts follows in the next section. 

3.3 Summary of the Six Gust Front Cases 

Following Colmer (1971), the parameters AV (magnitude of wind 
increase) and At (length of time interval over which the wind increases) 
are defined as the gust size and gust length respectively. If the speed 
of the advancing gust front, C, is constant, the spatial distance from 
the wind minimum to the wind maximum may be found by multiplying At by C . 
Several of the important parameters which characterize the six gust front 
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Figure 6. Time-height cross section of the hand- smoothed wind speeds for the Oklahoma 
gust front of April 26, 1969. The isotachs in this figure are expressed 
in mps. 







cases are presented in Table 1. These parameters include: 1) maximum 

wind speed; 2) gust size at 18 m; 3) gust size at the height of the 
maximum wind; 4) gust length at 18 m; 5) gust length at the height of 
the maximum wind; 6) the slope of the peak gust axis. 

In spite of the differences among the six gust front cases, several 
similarities may be noted. In five out of six cases, the gust size in- 
creased with height and the gust length decreased with height. These 
results agree with the average properties found by Colmer (1971). 

In all six cases, the axis of the maximum wind was nearly vertical. 

Both slight upstream and downstream slopes existed, but no preferred direction 
was evident. In all cases, the vertical wind shear increased over the gust 
length and reached a maximum at the time of maximum wind. 

The variation from case to case in slope of the peak gust axis as 
well as the gust front itself and the variation in gust length with 
height are probably results of surface friction. Colmer suggested 
that a lag in the gust front near the ground occurs intermittently, 
so that the denser cold air aloft spreads out ahead of the surface 
portion of the gust front. Turbulent mixing tends to destroy the 
overhanging bulge as the dense air sinks toward the surface. The lag 
then reforms through the effect of surface drag. These two non-steady 
processes in the immediate vicinity of the gust front may lead to small 
slopes in either the upstream or downstream direction at any time, and 
they are unpredictable for a particular storm at a particular location. 
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4.0 MODELING THE THUNDERSTORM GUST FRONT 


4.1 Scaled Time-Height Cross Sections 

Differences in magnitudes of the gust size, gust length, and height 
at which the peak gust occurs make it difficult to draw general conclusions 
about the vertical structure of the wind speed associated with the six gust 
fronts. Scaled time-height cross sections were constructed in order to 
test the hypothesis that the gust fronts differ only in magnitude, but 
not in basic structure. Empirical scaling parameters for time, height, 
and velocity were used to construct non-dimensional time-height cross 
sections of velocity. For the time scaling parameter, an obvious choice 
is gust length At. The height scaling parameter was chosen to be the 
height at which the maximum wind occurred, AZ. The velocity scaling 
parameter was chosen to be the gust size, AV. Of the three scaling 
parameters, AZ was the most difficult to define. As indicated in the un- 
sealed cross sections (Figures 2a, 3a, 4a, 5, 6, and 7), the maximum wind 
is elongated in the vertical, so that a precise definition of AZ is 
difficult. Also, in one Cape Kennedy case, the height of the maximum 
occurs above the top level of the tower; the height in this case could 
be estimated only approximately. 

Figures 8, 9, and 10 show the scaled cross sections for the three 
Cape Kennedy cases. In contrast to the unsealed cross sections, the 
scaled gust front profiles are quite similar. The maxima of scaled 
velocities are all about 1.7 and occur at Z/AZ and t/At of 1.0. These 
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Figure 10. Scaled time-height cross section of wind speed 
for the Florida gust front of June 25, 1968. 
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results indicate that if the scaling parameters AV, At, and AZ are known, 
at least part of the wind structure of the gust front could be deter- 
mined. The difficulty lies in predicting the scaling parameters, especially 
AZ. Considerable effort has been directed toward the prediction of 
peak surface wind gusts (which should be closely related to AV) . The 
parameter At might be related to parameters such as cell size and speed 
of movement, which could be determined from radar. However, there is 
no obvious basis on which to choose AZ. For this reason, the effort 
to find a height scaling parameter was abandoned and a simple hypothesis 
to explain the wind structure was tested. This hypothesis was simply 
that the vertical variation of the wind at the time of maximum wind 
speed might follow the logarithmic wind law in the lower portion of the 
gust front, because logarithmic profiles are known to occur with strong 
winds (Blackadar et al. , 1972). 


4.2 Test of the Logarithmic Wind Law Hypothesis 

The logarithmic wind law is derived with the assumption that Richardson 
number, is small, where 


Ri 


3 
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80 

3z 


c — ) 2 


(— ) 2 


(4.1) 


Here u and v are horizontal velocity components in arbitrary Cartesian 
directions. The Richardson numbers associated with the peak gust are 
undoubtedly small, regardless of the static stability, because of the 
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great wind shear. The logarithmic wind law states that the ratio of 
wind speeds (V^ and V^) at two different heights (z^ and z^) is a function 
only of the roughness length, z^. 


V 2 ln( 2 2 /z Q ) 


(4.2) 


The roughness length depends upon the upwind terrain. Therefore, at any 
one location, z^ is a function of wind direction. Recent results 
(Blackardar et al. , 1972) indicate that there is no systematic departure 
from the logarithmic wind law tjp to a height of about 150 m, which is 
higher than previously thought . 

The first step in testing the validity of the logarithmic wind law 

was the construction of time-height cross sections of the ratio 

(V/V- 0 ) of the wind speed to the speed at 18 m. If the logarithmic 
lo 

wind law is valid, then this ratio should be a function only of the 
roughness length, which may depend on wind direction, and also is 
different at Cape Kennedy and Oklahoma. In the construction of these 
cross sections, the time axis was scaled by At with the origin chosen 
as the time of the wind minimum prior to the gust front. Thus, the peak 
gust always occurs at t/At = 1. The 18 m wind speeds for Oklahoma were 
obtained by linear interpolation between the 7 m and 44.5 ih levels. 

The resulting isopleth analyses are shown in figures 11-16. 
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Figure 13. Ratio of wind speed to speed at 18 m vs. 

height and scaled time for the Florida 
gust front of June 25, 1968. 
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Figure 15. Ratio of wind speed to speed at 18 
height and scaled time for the Okl 
gust front of April 26, 1969. 




Figure 16. Ratio of wind speed to speed at 18 m vs, 
height and scaled time for the Oklahoma 
gust front of May 31, 1969. 
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The scaled velocities associated with the six scaled cross 


sections show several important similarities . The vertical shear is 

concentrated below 100 m, and the shear is relatively independent of 

time. Above 100 m, the structure is more complicated. On the average, 

the vertical shear is less. The several maxima and minima of shear 

above 100 m is most of the sections are not associated with the maxima 

and minima in acutal wind speeds. 

The scaled speeds in the three Oklahoma cases (Figures 14-16) 

show a greater range in the shear than the Florida cases. However, the 

range in the lower 150 m is about the same. Throughout the entire time 

interval, the shear in the lower levels is remarkably constant. 

The nearly steady-state ratio of V/V- 0 in the lower 100 m of each 

lo 

gust front supports the tenative hypothesis that the logarithmic wind 
law is valid in this region. The next step in testing this hypothesis 
was to determine whether the differences in the ratios among the six 
cases could be explained on the basis of different roughness lengths. 

4.3 Comparison of Roughness Lengths Under Thunderstorm Conditions 
With Those Under Non-thunderstorm Conditions 
In previous studies, Blackadar et al. (1972) and Sanders and 
Weber (1970) have empirically determined roughness lengths for the 
terrain around the Cape Kennedy and WKY Towers. Both studies were made 
under strong wind, but non- thunderstorm, conditions. Plots of 
the average z^, as a function of azimuth, from these studies are 
shown in Figures 17 and 18 for Cape Kennedy and Oklahoma. 
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Figure 17. Roughness length, as a function of direction at 

the Cape Kennedy tower (Blackadar et al., 1972). 
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the WKY Oklahoma tower. 


At Cape Kennedy (Figure 17), varies from about 0-3 m for 
northeasterly winds to about 0.8 m for southwesterly winds. In 
Oklahoma, the roughness lengths are much smaller than the Cape Kennedy 
values, ranging from about 0.01 m for easterly winds to about 0.14 m 
for southeasterly winds. 

To compare the values of that are computed under thunderstorm 
conditions to the previously determined values, ratios of the wind speed 
at 90 m to the speed at 18 m were averaged over the gust front for each 
case. Roughness lengths were then computed from 


lnz 


0 




lnz^ 


(4.3) 


The average ratios, computed roughness lengths, and average wind 
direction for each gust front, as well as previously determined rough- 
ness lengths, are listed in Table 2. These values of z^ are also shown 
in Figures 17 and 18 as a function of wind direction- The agreement between 
the values of z^ computed under thunderstorm conditions with the values 
determined under non-thunderstorm conditions is reasonably good in five 
out of six cases. This agreement implies that the logarithmic wind law 
is valid up to about 100 m during the gust front passage and that the 
differences between the speed ratios in the six gust front cases might 
be explainable on the basis of different roughness lengths associated 
with varying terrain and wind directions. 

The one Oklahoma case which has a roughness length of .290 m, 
compared to Sanders and Weber’s value of .022 m for the same wind direction, 
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suggests that more data must be analyzed in order to determine the 
best estimate of roughness length for a particular wind direction at 
the WKY Tower. The .022 m roughness length is based on a sample of 
only two cases, and when combined with a value of .290 m found in this 
study, indicates that the actual roughness length is probably between 
.22 m and .290 m. 

The directional variation of roughness lengths at a site is 
usually unknown, but by studying the surrounding terrain, a good 
qualitative estimate of the average roughness length may be made. Such 
an estimate could then be used to evaluate the ratios of wind speeds at 
high and low level. Therefore, an average ratio of upper level to lower 
level winds was determined separately for the three Cape Kennedy and 
the three Oklahoma cases. Each average ratio corresponds to an average 
roughness length for each location. 

Figures 19-21 indicate the observed time variation of the 90 m wind 

and the time variation of the 30 m wind profile multiplied by the 

average ratio equal to 1.25 for the three Cape Kennedy cases. 

Even though z does vary with direction, the use of the average ratio 
0 

to compute the upper level speed (implying a single average of 0.4 m) 
succeeds reasonably well. 

A similar comparison of upper level wind (90 m) vs. the lower 
level wind (7 m) multiplied by the average ratio of the winds at the 
two levels is shown in Figures 22-24 for the Oklahoma cases. 

This ratio was taken as 1.50, implying an average 2, of .04 m. 

Again, the agreement is quite good. These results indicate that if 
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Figure 19. Comparison of observed 90 m wind with the 30 m wind multiplied by the 
average ratio, v/v for the Florida gust front of July 27, 1967. 


FLORIDA 
April 29, 1968 


U 

O 

IO 


T° 



o o o o 

IO CM — 


(sdui) Q33dS 


43 


Figure 20. Comparison of observed 90 m wind with the 30 m wind multiplied by the 
average ratio, V qn /V n , for the Florida gust front of April 29, 1968. 



(sdui)Q33dS 


44 


Figure 21. Comparison of observed 90 m wind with the 30 m wind multiplied by the 
average ratio, v /V n , for the Florida gust front of June 25, 1968. 
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Figure 22. Comparison of observed 90 m wind with the 7 m wind multiplied by the 
average ratrio, V qr) /V , for the Oklahoma gust front of April >16 , 1969 
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Figure 23. Comparison of observed 90 m wind with the 7 m wind multiplied by the 

average ratio, V /V , for the Oklahoma gust front of April 26, 1969. 
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Figure 24. Comparison of observed 90 m wind with the 7 m wind multiplied by the 
average ratio, V on /V 7 , for the Oklahoma gust front of May 31, 1969. 


the time dependent behavior of the wind at one lower level is known, 
the variation of the wind at the other levels (up to 100 m) may be 
closely estimated. 

4.4 Implications of the Logarithmic Wind Law for the Prediction of 
the Gust Front Wind Speeds 

Because of the validity of the logarithmic wind law up to about 100 m 
during the passage of a thunderstorm gust front, the winds up to 100 m are 
related to a lower level wind in a simple way. Although the wind speed 
structure above 100 m is quite variable from case to case, the shear is 
much less than that below 100 m. Thus, as a first approximation, the 
gust front may be modeled by the logarithmic wind law below 100 m and zero 
shear above 100 m. Therefore, if the local roughness length and the time 
dependent wind speed at any one level are known, an estimate of the 
vertical cross section of wind speed can be made. 

The key parameters in specifying the wind speed behavior at a 
single level are the gust length, At, and the gust size, AV. If these 
parameters are known, as estimate of the time dependent wind speed 
profile may be made by assuming some functional form of the wind 
increase (for example, linear)over the interval At. 

The implication of the above results is that a crude estimate of the 
time-dependent, vertical wind speed structure during the passage of a 
gust front can be made from knowledge of three parameters: 1) the average 

roughness length, z^, which must be estimated for each location; 2) the 
gust size, AV, and; 3) the gust length, At. 
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The next question to be investigated is whether or not the latter 
two parameters can be statistically inferred from measurable synoptic 
scale and radar variables to enable short range predictions of the para- 


meters to be made. 



5.0 STATISTICAL RELATIONSHIPS OF GUST FRONT PARAMETERS 


TO RAWINSONDE AND RADAR VARIABLES 


5.1 Choice of Independent Variables 

The results of the previous sections have indicated that a crude 
model of the time-dependent vertical variation of the wind speed during 
the passage of a thunderstorm gust front may be constructed from knowledge 
of the roughness lenght 2^, the gust length, At, and peak gust, 

It is therefore important to determine whether or not these parameters are 
statistically related to dynamic and thermodynamic variables that could 
be obtained from nearby rawinsonde and/or radar observations. For this 
purpose, rawinsonde, radar, and surface observations for 81 thunderstorm 
cases were obtained from Tampa, Florida. Tampa was chosen because of 
the availability of nearby radar and rawinsonde data and also because 
thunderstorms there should be representative of subtropical thunderstorms. 
The variables chosen to serve as predictors for the parameters V 

max 

and At included: 1) the dry instability index, I; 2) maximum cloud 

height at the time of the peak gust; 3) environmental wind speed at 
the maximum cloud height; 4) maximum wind speed between 1000 and 500 
millibars; 5) the average surface wind speed prior ot the gust; 6) cell 
speed at the time of the maximum gust. 

The dry instability index, I, is a measure of the static stability 
of the atmosphere. Miller (1967) has previously shown that I is a useful 
predictor of peak gusts for Florida thunderstorms. The maximum cloud 
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height is indicative of the intensity of convection and should therefore 

be positively correlated with V . Because the downdrafts in cumulonimbus 
r max 

clouds transport upper level momentum to the surface, V should be 

max 

positively correlated with the wind speed at some upper level. Two measures 
of the upper level wind were tested as predictors; one, the wind at the 
cloud top level, and the other, the maximum wind speed below 500 millibars. 

The average wind speed at the surface prior to the gust was chosen 
as an independent variable for the regression analysis because it is 
plausible that this ambient velocity is an additive component of the 
gust velocity. 

5.2 The Variation of At Based on a Simple Kinematic Model 

The physical variables that determine At are more closely dependent 
upon the structure of the gust front rather than on larger scale variables. 
A schematic drawing of the thunderstorm and its associated gust front 
(Figure 25) illustrates some of the variables upon which At depends. 

We assume, based on observations, that the gust front precedes the pre- 
cipitation echo by a distance, Ar, that depends on the age of the storm. 

In general, the gust front is propagating outwardly with a radial 
velocity V , relative to the precipitation cell, which may also be 
translating with a velcoity C. If we assume that the gust length is 
proportional to the time between the gust front passage and the arrival of 
the precipitation echo, we may relate At to: 1) the perpendicular 

distance, from the station to the track of the storm center; 2) the 
radius, r, of the precipitation echo; 3) the speed of translation of the 
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Figure 25. Schematic diagram of gust front. 


precipitation cell; 4) the radial velocity of the gust front relative 
to the precipitation cell; and 5) the age, t, of the gust front, where 
t=0 is defined as the time when Ar=0. This relationship is 


V t sin (90 + J [arcsin (-) - arcsin (— — ) ]) 

At = ( — —) { £ - v rt+r 

^ C ; 1 r . ,1, . , a 


sin (90 - — [arcsin (— ) + arcsin ( TT _ , 
l r V-rt+r 


) ]) 


(5.1) 


Figure 26 shows the dependence of At on the age of the cell and the 
radial velocity of the gust front. Each curve was determined for a 
station 800 m from the center of the storm track for a precipitation cell 
of 1000 m radius translating at 5 mps. The three curves corresponding to 
different radial velocities indicate that At is very sensitive to the 
radial velocity of the gust front over the gust front lifetime, so that, 
a slight increase in radial velocity produces a large increase in At with 
time. 

The variation of At with precipitation cell radius and station dis- 
tance from the center of the storm track is illustrated in Figure 27. 

In this calculation the precipitation cell speed, radial velocity of gust 
front and gust front age were defined by C=5 mps, V^=2 mps, and t=10 
minutes. As shown in Figure 27, At does not vary with r at £=0. Further- 
more, At is relatively insensitive to the distance of the station from 
the center of the storm track until that distance exceeds about two-thirds 
of the cell radius. In each case, At increases very rapidly as Z approaches 
the cell radius r. Finally, for a given 5,^0, the storms with the smaller 
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Figure 26. At vs. gust front age (Time). 
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Figure 27. At vs. station distance from center of storm 
track (&) . 


radii are associated with the largest values of At. 

In summary, a simple kinematic model yields reasonable values of 
gust length At as a function of the age of the storm, the speed of the 
storm movement, outward propagation velocity of gust front with respect 
to the thunderstorm echo, cell radius, and distance of the station from 
the track of the storm. For a typical range of these parameters, At 
varies from zero to ten minutes. 

5.3 Development of Regression Equations 

The regression analyses were performed on an IBM 360/67 computer. 

Means and standard deviations of each variable and correlation coefficients 
for each possible pair of variables were computed. These means, standard 
deviations, and correlation coefficients are listed in Table 3. 

The first column of Table 3 lists the independent and dependent 
variables. The second and third columns give the respective means and 
stantard deviations for each variable. The last six columns, when 
matched with the eight variables in column one, show the correlations 
between all possible pairs of variables. 

As expected, the cell speed is positively correlated (r=0.37) 
with the maximum wind in the 1000-500 millibar layer. A relatively strong 
correlation (0.32) also exists between the instability index and the 
average surface wind speed prior to the gust. The higher average wind 
speeds are associated with the more unstable conditions. The physical basis 
for this correlation is the increased vertical exchange of momentum under 
unstable conditions. Under unstable conditions, therefore, high momentum 
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from aloft is carried downward and low momentum near the surface is 
transported upward. This reasoning also explains the significant 
negative correlation of “0.27 between the instability index and the wind 
speed at cloud tops. 

The second part of the analysis was a step-up multiple linear 
regression procedure, in which number of multiple linear regression equations 
were computed. These equations are obtained by adding or deleting one variable 
at a time, giving regression equations at intermediate steps 


o o 
II 

+ 

b J x i 



rH O 

rO 

II 

+ 

b i x i 

+ 

b 2 X 2 

CN O 
II 

+ 

b J x i 

+ 

b 2 x 2 + b 3 x 3 


/s i n~l , - n-1 , , n-1 , 

y n ■ b 0 + b l X 1 + b 2 X 2 + 


. , n-1 

+ b x 
n n 


(5.2) 


where the y_^ are the best estimates in a least squares sense of the 
independent variable (the jjredictand) and the x^ are the various predictors. 
The coefficients, b^, in equation (5.2) are determined by minimizing 
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E = E (y - y . )‘ 
1 


(5.3) 
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Table 3. Means and standard deviations for each variable and correlation coefficient for 

each possible pair fo variables for the 81 thunderstorm cases from Tampa, Florida. 
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the sums of squares of the deviations of the observed value, from 
the predicted value y^. In determining the intermediate regression 
equations, Variables were added in the order which yielded the greatest 
reduction in the error E. If the variance reduction associated with 
a given independent variable is insignificant at a specified level 
(2.5 F level), the variable is not included in the regression. There- 
fore, the final equation contains only those variables that explain important 
portions of the variance. A more detailed discussion of the step-up 
regression procedure is given by Efroymson (1967). 

5.4 Discussion of Regression Equation 

The results from the 81 Tampa thunderstorm cases showed that no 
independent variables were significantly correlated with the para- 
meter At and no regression equation was formed. This result indicates 
that At cannot be estimated from the larger scale rawinsonde and radar 
parameters studied here. As discussed in section 5.2, At is a strong 
function of variations in cell radius, cell age, motion of the cell, and 
position of the station with respect to the center of the storm track, 

Since these variables were unknown, it is not suprising that At was un- 
correlated with any of the independent variables. 

In the analysis for the peak wind gust, two independent variables, 
the surface wind speed prior to the gust and the dry instability index, 
were correlated relatively strong with the peak gust. Figure 28 shows a 
scatter diagram of the gust speed vs. the average surface wind speed. 
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Figure 28, Scatter diagram with line of best fit of the 
speed maximum vs. the average speed before 
the gust front. 
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The line of best fit is also shown. Figure 29 shows the gust speed as a 
function of the instability index, with the line of best fit also shown. 

Both figures exhibit considerable scatter. 

The regression equation relating the peak gust (in mps) to the 
average surface wind and the instability index is 

V = 3.2 + 1.02 V + .158 I (5.4) 

max s 

where V is expressed in mps and I is expressed in degrees Celsius, 
s 

However, these two independent variables explained only thirty-four 
percent of the variance of peak gusts, and therefore are of only limited 
use in prediction. 

The small percentage of explained variance for is somewhat 

disappointing, but not unexpected in view of the complexity of the unsteady, 
three-dimensional thunderstorm circulation. There are several possible 
explanations why so little variance was explained. First, there may 
be additional important independent variables that were not tested. 

The possibility of finding additional significant variables is considered 
unlikely because environmental variables that are known to be important 
in thunderstorm development (static stability and vertical wind shear) 
have already been included. 

A second and more plausible explanation is that the rawinsonde 
parameters were not representative of the environment in the immediate 
vicinity of the thunderstorm cell. Many of the rawinsonde data were six 
hours earlier or later that the thunderstorm, during which time the 
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Figure 29, Scatter diagram with line of best fit of the 
speed maximum vs. the instability index. 
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environmental variables could have changed appreciably. It is probable 
that simultaneous soundings in the immediate vicinity of the thunderstorms 
would provide further reduction of the scatter. 

A third possibility is that dependent variable V^ x is a point 
estimate and may not be representative of the average peak gust associated 
with the thunderstorm. In the three-dimensional gust front structure, 
wide variation with respect to age of the gust front, and position with 
respect to the thunderstorm center may be expected. Thus, to expect any 
real reduction in scatter, one would have to stratify the data according 
to these two parameters. 
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